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1  5th  Year  Executive  Summary 

This  report  summarizes  the  fifth  year  of  the  MURI  in  RF  Photonic  Systems  which  assem¬ 
bles  the  University  of  Colorado,  the  University  of  California-Davis,  Montana  State  University, 
and  George  Mason  University  into  a  systems  oriented  research  team  investigating  the  ap¬ 
plication  of  photonic  techniques  to  the  control  and  processing  of  RF  phased  arrays.  Under 
this  MURI  funding,  the  prime  contract  is  administered  through  the  University  of  Colorado 
under  the  direction  of  Prof.  Kelvin  Wagner. 

Photonic  techniques  have  emerged  as  the  preferred  approach  to  RF  communication  and 
remoting  tasks  throughout  the  microwave/millimeter-wave  bands,  when  size,  weight,  power, 
and  dispersion  free  bandwidth  are  key  factors.  In  addition,  applications  such  as  controlling 
time  delay,  spectral  filtering,  mixing,  frequency  generation,  antenna  array  control,  beam 
forming,  wide-band  signal  processing,  and  target  recognition  are  being  implemented  opti¬ 
cally.  Much  of  this  development  has  been  driven  by  the  evolution  of  high  frequency  photonic 
devices  such  as  100  GHz  detectors,  40  GHz  modulators,  20  GHz  lasers,  and  3  THz  fiber  optic 
transmission  and  delay  lines.  Often  these  technologies  have  been  used  as  direct  replacement 
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of  individual  RF  components  with  their  optical  counterparts,  relying  on  the  smaller  size 
and  lighter  weight  of  the  optical  components  and  waveguides  to  yield  a  resulting  system 
advantage.  More  sophisticated  applications  utilize  the  incredibly  large  bandwidth  of  opti¬ 
cal  transmission  to  simultaneously  wavelength-multiplex  a  large  number  of  RF  signals  for 
efficient  parallel  transmission  or  to  control  the  element  time  delays  through  dispersion.  The 
approach  of  this  MURI  team  takes  RF  photonics  to  an  even  higher  level  of  sophistication,  in 
which  the  massive  parallelism  possible  using  dynamic  volume  holography  in  photorefractive 
and  photon-echo  materials  enables  the  implementation  of  optimal  adaptive  and  nonlinear 
algorithms  for  array  processing  in  the  optical  domain. 

The  system  research  groups  comprising  this  MURI  team  have  pursued  several  differ¬ 
ent  avenues  towards  one  common  goal  of  inserting  photonics  into  RF  array  systems.  We 
have  developed  solutions  that  build  on  the  strengths  of  each  approach,  and  identified  the 
optimal  combination  of  devices,  materials,  algorithms,  and  systems.  The  BEAMTAP  ap¬ 
proach  (Broadband  and  Efficient  Adaptive  Method  for  True-time-delay  Array  Processing) 
has  been  investigated  theoretically  and  implemented  experimentally  to  demonstrate  efficient 
adaptive  beam  forming  using  acousto-optic  deflectors,  traveling-fringes  detectors,  and  pho¬ 
torefractive  adaptive  weights.  The  novel  physics  of  optical  coherent  transients  (OCT)  have 
been  developed  for  applications  as  programmable  RF  time  delays  as  well  as  array  signal 
processors,  since  these  nonlinear  materials  can  directly  produce  the  necessary  delayed  RF 
signals  as  photon  echoes.  Both  of  these  approaches  utilize  coherent  beam  forming  in  the 
optical  domain  which  can  improve  the  noise  figure  of  optically  remoted  arrays  and  opti¬ 
cal  novelty  filtering  developed  at  CU  can  enable  robust  RF  signal  combination  even  in  the 
presence  of  phase  errors  and  drifts.  Optically  controlled  quasi-optical  active  antenna  arrays 
provide  a  near  term  test  vehicle  for  optically  controlled  arrays,  allowing  rapid  switching 
between  polarizations,  transmit /receive  (T/R)  mode,  and  frequency  response,  and  in  addi¬ 
tion,  demonstrate  the  capability  of  producing  multiple  simultaneous  true-time-delay  (TTD) 
beams  as  a  front  end  for  beam-space  adaptive  optical  processors.  Optical  nonlinearities  have 
been  explored  as  an  approach  to  dynamically  solving  the  problem  of  signal  extraction  from 
small  RF  antenna  arrays.  Practical  polymeric  in-line  fiber  (PILF)  modulators  with  high¬ 
speed  coplanar  waveguide  (CPW)  electrode  structures  have  been  developed  as  the  front  end 
RF-to-optical  transducers  for  the  antenna  arrays.  All  of  these  efforts  have  been  tied  to¬ 
gether  by  investigations  of  spatio-temporal  signal  processing  algorithms  compatible  with  RF 
photonic  technology,  which  allows  us  to  evaluate  and  compare  the  various  optical  systems 
techniques  being  investigated.  As  an  example  system  developed  as  a  unique  collaboration 
between  MURI  team  members,  an  RF  antenna  array  has  been  coupled  into  a  photorefractive 
optical  system  for  extracting  the  principle  component  in  the  RF  signal  environment,  and  the 
resulting  system  has  been  packaged  as  a  portable  suitcase  demonstrator. 
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2  Education,  Presentations,  Outreach,  and  Manage¬ 
ment 

2.1  MURI  Team 

This  MURI  effort  has  directly  or  indirectly  funded  the  6  principal  investigators,  5  post¬ 
doctoral  or  visiting  research  scientists,  17  graduate  students,  11  undergraduate  students,  2 
research  assistants,  and  fractions  of  a  financial  manager  and  secretary. 

2.1.1  Principal  Investigators 

The  original  MURI  team  of  6  Pis  was  joined  by  Andre  Knoesen  from  UC  Davis  in 
year  2  year,  replacing  Ted  Weverka.  Prof.  Andre  Knoesen  worked  on  PILF  modulators 
on  a  separately  funded  subcontract  which  ends  this  year.  Weverka’s  project  on  coherent 
optical  beam  combination  for  phased  array  processing  is  being  continued  at  the  University 
of  Colorado  as  the  required  front  end  to  the  BEAMTAP  processor. 


Table  1:  Principal  investigators  involved  in  the  RF  Photonic  Systems  program. 


Researcher 

Department 

Institution 

Prof.  Kelvin  Wagner 

Prof.  Zoya  Popovic 

Prof.  Dana  Anderson 

Prof.  William  Randall  Babbitt 
Dean  Lloyd  Griffiths 

Prof.  Andre  Knoesen 

ECE/OCS 

ECE 

JILA 

Physics 

ITE 

ECE 

CU-Boulder 

CU-Boulder 

CU-Boulder 

Montana  State  University 
George  Mason  University 
UC- Davis 

2.2  Students,  Post-Doctoral  Researchers,  and  other  Personnel 
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Table  2.  Undergraduate  students  who  have  been  involved  in  MURI  related  research  dur¬ 
ing  the  past  year.  Again,  those  who  worked  on  closely  related  projects  in  collaboration 
with  the  MURI  effort  are  indicated  by  *. 


Undergrad 

Status 

Position 

Advisor 

Dept. 

Inst. 

Eric  Hoyt 

graduated  BS/MS  12/02 

Undergraduate 

Wagner 

OCS 

CU 

John  Jost 

Assistant 

Undergraduate 

Anderson 

Phys 

CU 

Joe  Fischer* 

graduated  6/02 

Undergrad 

Babbitt 

Physics 

MSU 

Zeb  Barber 

graduated  6/03 

Undergrad 

Babbitt 

Physics 

MSU 

Jesse  Law 

continuing 

Undergrad 

Babbitt 

Physics 

MSU 

Table  3:  MURI  personnel  indicating  those  funded  by  this  program.  Those  who  worked 
on  closely  related  projects  in  collaboration  with  the  MURI  effort  are  indicated  by  *. 


Researcher 

Status 

Position 

Advisor 

Dept. 

Inst. 

Andrew  Kiruluta 

6/02-11/02 

Post-doc 

Wagner 

OCS 

CU 

Steve  Jia 

continuing 

Grad  student 

Babbitt 

Physics 

MSU 

Mingzhen  Tian 

continuing 

Research  Scientist 

Babbitt 

Physics 

MSU 

Diego  Yankelevich 

3/99- present 

Adjunct  Prof. 

Knoesen 

EE 

Davis 

Joe  Shamir 

9/02-10/02 

Visiting  MURI  fellow 

Wagner 

OCS 

CU 

Edilene  Fotheringham 

continuing 

Grad  Student 

Anderson 

JILA 

CU 

Greg  Kriehn 

DoD  fellowship 

Grad  Student 

Wagner 

OCS 

CU 

Friso  Schlottau 

continuing 

Grad  Student 

Wagner 

OCS 

CU 

Paul  Smith* 

continuing 

Grad  Student 

Popovic 

ECE 

CU 

Jan  Peeters  Weem 

PhD  finished  6/02 

Grad  Student 

Popovic 

ECE 

CU 

Jacques  Hong  Loui 

continuing 

Grad  Student 

Popovic 

ECE 

CU 

Bob  Peters* 

MS  finished  8/01 

Grad  Student 

Babbitt 

Physics 

MSU 

Randy  Reibel 

Graduated  8/02 

Grad  Student 

Babbitt 

Physics 

MSU 

Hongyan  Li 

continuing 

Grad  Student 

Babbitt 

Physics 

MSU 

John  Campbell 

continuing 

Grad  student 

Knoesen 

EE 

Davis 

Liu  Ming  Wu 

continuing 

Grad  student 

Knoesen 

EE 

Davis 

Leslie  Czaia 

continuing 

Prof.  Res.  Asst. 

Anderson 

JILA 

CU 

John  Getty 

EE  engineer 

Babbit 

Physics 

MSU 

Norm  Williams 

Machinist 

Babbit 

Physics 

MSU 
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2.2.1  ^  Students  Graduated  from  the  MURI  program 

1.  Jan  Peeters  Weem,  June  2002,  CU,  ECE,  PhD,  Popovic.  Thesis  title:  “Broadband 
antenna  arrays  and  noise  coupling  for  radio  astronomy.” 

2.  Randy  Reibel,  August  2002,  MSU,  Physics,  PhD,  Babbit. 

3.  Edeline  Fotheringham,  May  2003,  CU,  ECE,  PhD,  Anderson.  Thesis  title:  “Microwave 
signal  processing  with  photorefractive  dynamic  holography.” 

4.  Greg  Kriehn,  May  2003,  CU,  ECE,  PhD,  Wagner.  Thesis  title:  “Coherent  optical 
signal  processing  for  broadband  adaptive  phased-array  antennas  using  the  BEAMTAP 
algorithm.” 

2.3  Meetings 

The  Boulder  team  met  weekly  to  exchange  information,  make  presentations,  as  well  as 
to  discuss  technical  and  management  details.  The  entire  MURI  team  gathers  once  or  twice 
a  year  for  intellectual  exchange,  discussions,  and  review  regarding  scientific  progress  during 
the  year.  One  of  these  meetings  is  primarily  an  internal  MURI  meeting,  and  the  other  is  a 
formal- review  for  programmatic  purposes.  Our  third  annual  review  was  held  in  Boulder  on 
October  12  2000  conjunction  with  the  UCLA  MURI  review  on  October  13.  Plans  for  the 
4th  annual  review  to  be  held  in  Washington  DC  on  Nov  30,  2001  were  put  on  hold  when  the 
Program  Manager,  Dr.  Steve  Pappert  of  SPAWAR,  took  a  leave  from  the  government. 

2.3.1  5th  year  MURI  team  Interactions 

1.  Prof.  Popovic’s  group  has  collaborated  mainly  with  Prof.  Anderson’s  group  at  CU, 
and  somewhat  with  Prof.  Griffiths  at  MSU. 

2.  Prof.  Popovic  visited  Prof.  Knoessen’s  group  at  UC  Davis  and  gave  a  graduate  semi¬ 
nar,  October  2002. 

3.  Prof.  Babbitt  visited  UC-Boulder  in  November  of  2002. 

2.4  Additional  Collaborations  and  Related  Research 

The  MURI  has  enabled  a  wide  variety  of  collaborations  among  the  Pis  as  well  as  with 
other  groups.  In  this  section  some  notable  collaborations  are  discussed  and  outgrowths  from 
the  MURI  into  new  areas  are  mentioned. 
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1.  Prof.  Joseph  Shamir  of  the  Technion  in  Israel  was  the  MURI  fellow  and  visited  Prof. 
Wagner’s  group  between  Aug.  and  Sept,  of  2001  and  between  Sept,  and  Oct.  of  2002. 

2.  Prof.  Vitaly  Voloshinov  of  Moscow  State  University  in  Russia  was  the  MURI  fellowship 
visitor  to  Prof.  Wagner’s  group  to  work  on  acousto-optic  devices,  Nov  11-21,  2001. 

3.  Prof.  Babbit  explored  Raman  amplification  at  793nm  for  OCTs  with  Directed  Energy 
Systems  of  Colorado  Springs  in  the  Fall  of  2002. 

4.  Two  graduate  students  under  the  supervision  of  Prof.  Popovic  have  been  for  3  months 
at  the  NFRA  facility  at  Dwingeloo  (Jan  Peeters  Weem  and  Eric  Bryerton).  Several 
papers  have  resulted  from  this  collaboration. 

5.  Technishe  Universitat  Miinchen  (Germany)  and  Univ.  of  Karlsruhe  (Germany):  Two 
German  students,  Holger  Matern  and  Johannes  Russer  were  sent  (and  supported  by 
their  home  country)  to  work  in  Prof.  Popovic’s  group  on  their  MS  theses  in  the  area 
of  RF  photonics.  Holger  developed  a  5-channel  electrooptic  modulator  and  Johannes 
is  finishing  a  fast  optically-controlled  microwave  phase  shifter. 

2.4.1  Other  related  grants 

Profs.  Popovic  and  Wagner  submitted  new  proposals  on  optical  processing  of  RF 
signals  from  broadband  RF  antenna  arrays.  Any  subsequent  grant  will  be  a  direct 
result  of  the  Co-PI’s  joint  work  under  this  MURI. 

2.  Prof.  Popovic  and  Anderson  have  received  a  follow-on  grant  from  NSF  to  expand  the 
optical  processing  to  wireless  communications.  This  grant  is  now  ending,  and  they 
submitted  a  follow-on  proposal  to  the  NSF  ITP  program,  Dec.  2002.  The  proposal  is 
currently  under  review. 

3.  “High  Performance  Correlators  based  on  Spectral  Hole  Burning  Technology,”  NASA, 
PI:  Alan  Craig,  $5,096,000,  4/99-4/03.  PI  changed  to  W.  R.  Babbitt  in  12/2001. 

4.  Prof.  Babbit  received  an  award  from  Montana  Board  for  Research  and  Commercial¬ 
ization  of  Technology:  matching  grant  for  the  MURI  grant  $7,810,  5/01-5/02. 

5.  “Consortium  for  Stable  Laser  Applications,”  NSF  Partners  for  Innovation  (PFI),  PI: 
Alan  Craig,  co-PI  Jim  McMillan,  $598,553,  10/01-9/03,  changed  PI  to  W.  R.  Babbitt 
in  12/01. 

6.  “Applications  of  Optical  Coherent  Transient  Technology:  Pulse  Shaping,  Spectral  Fil¬ 
tering,  Arbitrary  Waveform  Generation,  and  RF  Beamforming,”  AFOSR/  DEPSCOR, 
PI:  W.  R.  Babbitt,  co-PI:  Mingzhen  Tian,  $310,067  plus  $155,034  match,  7/02-7/05. 
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7.  “Wide-Bandwidth  and  Multidimensional  Analog  Signal  Processing,”  DARPA,  Lead: 
W.  R.  Babbitt,  Co-PIs  Kelvin  Wagner,  Univ.  of  Colorado  and  Randy  Equall,  Scientific 
Materials  Corp.,  and  Kris  Merkel,  Krishna  Rupavatharam,  and  Rufus  Cone  at  MSU, 
$1.9M,  10/02-9/04. 

8.  “Spatial  Spectral  Coherent  Holographic  Integrating  Processor-S2-CHIP:  Device  char¬ 
acterization,  performance  evaluation,  and  optimization,”  Montana  Board  for  Research 
and  Commercialization  of  Technology,  PI:  W.  R.  Babbitt,  co-PIs:  Kris  Merkel  and 
Mark  Ivey, $149, 420,  11/02-7/03,  matched  by  NASA  grant. 

9.  Prof.  Wagner,  $60K  subcontract  to  University  Nebraska  Lincoln  for  ONR  proposal  on 
Random  Noise  Radar.  Title:  Acoustooptic  Delay  Line/Correlator  for  random  noise 
radar,  8/00  -  8/02. 

10.  $350K  Prof.  Wagner  PI,  Real-time  Multibeam  Imaging  and  High-Resolution  Spectral 
Analysis  for  Large,  Wideband  Antenna  Arrays  -  submitted  to  NRO  DII  (National 
Reconnaissance  Office  Directors  Innovation  Initiative).  $350,000  for  9  months.  Status: 
Awarded  on  Dec  21,  2001,  rescinded  award  due  to  NRO  budget  cut. 

11.  $20K  Prof.  Wagner  PI,  Joint  Optoelectronic  Project  Proposal  -  Double-clad  Thulium 
doped  fiber  for  use  in  a  fiber  amplifier  for  optical  coherent  transient  experiments.  An 
equipment  proposal  to  the  Joint  Optoelectronics  Project  (JOP).  Status  AWARDED: 
$20,000  in  custom  optical  fiber 

12.  $14.5K  Prof.  Wagner  PI,  NSF/CNRS  collaboration,  multidimensional  photon  echo 
optical  processing,  notice  of  award  Jan  3  2000,  in  effect  2000-2004. 

13.  Prof.  Wagner,  co-PI,  $80K  DEPSCOR  subcontract  from  MSU,  Optical  coherent- 
transient  time-delay  control  using  acoustooptic  distributed  local  oscillators,  $80K, 
plus  $40K  CU  matching  7/2002-12/2004. 

14.  Prof.  Wagner  PI,  Coherent  Optical  Array  Processing,  DARPA  Analog  Optical  Signal 
Processing  (AOSP)  program,  awarded  in  2002.  In  effect  2002-2006. 

2.5  Awards  and  Honors  for  the  Pis  and  students 

1.  Best  Poster  Award  granted  to  Friso  Schlottau,  Randy  Babbit,  Jean-Louis  Le  Gouet  and 
Kelvin  Wagner  for  the  CPIA  poster  entitled  ’’Spatial  spectral  holography  for  wideband 
spectral  analysis,”  Nov.  2002. 
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2.  MWP  2001  on  Jan.  2002.  Best  student  presentation  runner-up:  Friso  Schlottau  and 
Kelvin  Wagner,  RF  Photonics  for  Simultaneous  Multiple  TTD  Beamforming  for  2-D 
Antenna  Arrays,”  IEEE  Microwave  Photonics  Meeting,  Long  Beach,  Oct  2001  (delayed 
to  Jan.  2002). 

3.  Prof.  Kelvin  Wagner  was  elevated  to  the  status  of  fellow  of  the  OSA  (Optical  Society  of 
America)  for  Contributions  to  adaptive  spatio-temporal  optical  processing  systems.” 

4.  Zoya  Popovic  and  Dana  Anderson  both  received  the  prestiguous  Humboldt  Award  for 
Senior  US  Scientists  from  the  German  Alexander  von  Humboldt  Stiftung  (2001). 

5.  Zoya  Popovic  received  the  HP/ASEE  Terman  Award  for  combined  teaching  and  re¬ 
search  excellence,  October  2001. 

6.  Zoya  Popovic  became  a  Fellow  of  the  IEEE. 

7.  Zoya  Popovic  received  the  Outstanding  Speaker  Award  at  the  ONR  Workshop  on 
Challenges  for  Multifunctional  Digital  TVansmit  Arrays,  Marco  Island,  Nov.  2001. 

8.  Zoya  Popovic  received  the  CU-Subaru  Educator  Spotlight  Award,  2003. 

9.  Prof.  Babbit  was  awarded  the  Charles  and  Nora  L.  Wiley  Faculty  Award  for  Merito¬ 
rious  Research  and  Creativity,  2002. 

10.  Ph.D.  student  Randy  Reibel  received  MSU’s  Graduate  Achievement  Award,  2003. 

11.  Undergraduate  student  Zeb  Barber  was  nominated  a  Goldwater  Scholar  (2002)  and 
graduated  with  Highest  Honors  (6/2003). 

12.  Randy  Babbitt,  became  Director  of  Spectrum  Lab,  February  2002-present. 

13.  Committee  chairs  for  Randy  Babbitt:  Chair  of  program  subcommittee,  Holography, 
Wavemixing,  Photorefractives,  and  Storage,  Conference  on  Lasers  and  Electro-Optics, 
Long  Beach,  CA,  May  2002.  Chair  of  program  subcommittee,  Holography,  Wavemix¬ 
ing,  Photorefractives,  and  Storage,  Conference  on  Lasers  and  Electro-Optics,  Balti¬ 
more,  MD,  June  2003.  Co-Chair,  Spectral  Holeburning,  Single  Molecules,  and  Related 
Spectroscopies  2003,  Bozeman,  MT,  July  2003. 

2.6  MURI  Publications 

2.6.1  Papers  published  in  2002  (since  8/15/01) 

1.  Andrew  Kiruluta,  G.  S.  Pati,  Gregory  Kriehn,  Paulo  E.  X.  Silveira,  Anthony  W.  Sarto 
and  Kelvin  Wagner,  “Spatio-temporal  operator  formalism  for  holographic  recording 
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and  diffraction  in  a  photorefractive  based  true-time-delay  phased  array  processor,” 
accepted  for  publication  in  Applied  Optics,  Sept.  2003. 

2.  J.  Shamir  and  K.  Wagner,  “Generalized  Bragg  selectivity  in  volume  holography,”  Ap¬ 
plied  Optics  41(32),  pp.  6773-6785,  Nov.  2002. 

3.  Paulo  E.  X.  Silveira,  G.  S.  Pati  and  Kelvin  H.  Wagner,  “Optical  Finite  Impulse  Re¬ 
sponse  Neural  Networks,”  Applied  Optics,  vol.  41,  no.  20,  pp.  4162-4180,  2002. 

4.  Steve  Blair  and  Kelvin  Wagner,  Gated  Logic  with  Optical  Solitons,  in  “Collision-Based 
Computing”,  ed  by  Andrew  Adamatzky,  Springer- Verlag,  UK,  to  be  published  in  2002. 

5.  Paulo  E.  X.  Silveira,  G.  S.  Pati  and  Kelvin  H.  Wagner  “Optoelectronic  implementation 
of  a  256-channel  sonar  adaptive  array  processor,”  to  be  submitted  to  Applied  Optics, 
2003. 

6.  D.  Z.  Anderson,  V.  Damiao,  D.  Popovic,  Z.  Popovic,  S.  Romisch,  A.  Sullivan,  “-70dB 
optical  carrier  suppression  by  two-beam  coupling  in  photorefractive  media,”  Applied 
Physics  B,  72,  pp.  743-748,  2001 

7.  E.  Fotheringham,  S.  Romisch,  P.  Smith,  D.  Popovic,  D.  Anderson  and  Z.  Popovic, 
“A  lens  antenna  array  with  adaptive  optical  processing,”  IEEE  Trans.  Antennas  and 
Propagation,  Special  Issue  on  Wireless  Communications,  vol.  50,  no.  5,  pp.  607-617  , 
May  2002. 

8.  H.  Loui,  J.  Peeters  Weem,  Z.  Popovic,  “Dual-band  dual-polarized  nested  Vivaldi  slot 
array  with  multi-level  ground  plane,”  to  appear  in  IEEE  Trans,  on  Antennas  and 
Propagation,  2003. 

9.  Ja Peeters  Weem,  Z.  Popovic,  “Noise  coupling  in  active  phased  array  antennas,”  to  be 
submitted  to  the  IEEE  Trans,  on  Microwave  Theory  and  Techniques. 

10.  P.  Smith,  E.  Fotheringham,  Z.  Popovic,  D.  Anderson,  “An  optoelectronic  implementa¬ 
tion  of  independent  component  analysis,”  to  be  submitted  to  IEEE  Trans  on  Circuits 
and  Systems,  August  2003. 

11.  R.  Reibel,  Z.  Barber,  M.  Tian,  and  W.  R.  Babbitt,  “High  bandwidth  spectral  grating 
programmed  with  linear  frequency  chirps,”  J.  Lumin.  98,  355-365  (2002). 

12.  Z.  Barber,  M.  Tian,  R.  Reibel,  and  W.  R.  Babbitt,  “Optical  pulse  shaping  using  optical 
coherent  transients,”  Opt.  Exp.  10,  1145-1150  (2002) 

13.  R.  Reibel,  Z.  Barber,  M.  Tian,  W.R.  Babbitt,  Z.  Cole,  and  K.D.  Merkel,  “Amplification 
of  high  bandwidth  phase  modulated  signals  at  793nm,”  J.  Opt.  Soc.  Am.  B.  19,  pp. 
2315  (2002). 

14.  R.  Reibel,  Z.  Barber,  J.  Fischer,  M.  Tian,  and  W.  R.  Babbitt,  “High  Bandwidth  Linear 
Sideband  Chirped  Programming  for  Optical  Coherent  Transients,”  to  be  published  in 
2003. 
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15.  J.  Campbell,  A.  Knoesen  and  D.  R.  Yankelevich,  “Measurement  of  the  modulation 
efficiency  of  an  optical  phase-modulator  using  a  self-homodyne  receiver,”  IEEE  Phot. 
Tech.  Letters,  vol.  14,  no.  9,  pp.  1130-1132,  2002. 

2.6.2  Conference  Presentations  during  2002  (since  8/15/01) 

1.  G.  Kriehn  and  K.  Wagner,  “Experimental  Demonstration  of  Adaptive  Phased-Array 
Signal  Processing,”  (poster)  Colorado  Photonics  Industry  Association ,  November  2002. 

2.  G.  Kriehn,  F.  Schlottau,  and  K.  Wagner,  “Optically-Implemented  2-D  Beam  Steering 
and  Jammer  Nulling  using  BEAMTAP,”  Optical  Computing ,  April  2002. 

3.  G.  Kriehn,  F.  Schlottau,  G.  S.  Pati,  and  K.  Wagner,  “Demonstration  of  RF  Photonic 
Beam  Forming  using  the  BEAMTAP  Algorithm,”  General  Electric  Meeting ,  November 
2001. 

4.  G.  Kriehn  and  K.  Wagner,  “Experimental  Adaptive  Beam  Forming  with  Polarization 
Read-Write  Multiplexing  using  BEAMTAP,”  (poster)  In  The  International  Topical 
Meeting  on  Microwave  Photonics,  IEEE,  Jan.  2002. 

5.  G.  Kriehn,  “Optical  Phased  Array  Signal  Processing,”  (poster)  Colorado  Photonics 
Industry  Association,  November  2001. 

6.  G.  Kriehn  and  K.  Wagner,  “Experimental  Adaptive  Beam  Forming  with  Polarization 
Read-Write  Multiplexing  using  BEAMTAP,”  The  International  Topical  Meeting  on 
Microwave  Photonics  MWP  2001,  January,  2002. 

7.  Balakishore  Yellampalle,  Kelvin  Wagner,  and  Steve  Blair,  “Anti-guide  assisted  spatial 
soliton  logic  gate,”  OSA  topical  meeting  on  Non-linear  Optical  Guided  Waves,  Como, 
Italy,  Sept.  2002. 

8.  Kelvin  Wagner,  Friso  Schlottau,  and  Jaap  Bregman,  “Array  imaging  using  spatial- 
spectral  holography,”  IOG  Information  Optics  meeting,  Mannheim,  Germany,  Sept. 
2002. 

9.  Alexandre  R.  S.  Romariz  and  Kelvin  Wagner,  “Optoelectronic  implementation  of  a 
FitzHugh-Nagumo  neural  model,”  Neural  Information  Processing  System  (NIPS),  Van¬ 
couver,  CA,  Dec.  2002. 

10.  Max  Colice  and  Kelvin  Wagner,  “Phase-cohering  holography  for  coherent  analog  opti¬ 
cal  signal  processing,”  GOMAC  03,  Tampa  FL,  Apr.  2003. 

11.  K.  H.  Wagner,  M.  Colice,  G.  Kriehn,  F.  Schlottau,  and  R.  T.  Weverka,  “Photonics 
Multiple  Beam  Forming  for  Broadband  RF  Antenna  Arrays,”  GOMAC  03,  Tampa, 
FL,  Apr.  2003. 

12.  W.  R.  Babbit,  K.  D.  Merkel,  M.  Tian,  R.  Krishna  Mohan,  Z.  Cole,  Y.  Sun,  Rufus  R.  L. 
Cone,  K.  H.  Wagner,  and  R.  W.  Equall,  “Optical-Coherent-Transient  technologies  for 
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wide-bandwidth  and  multidimensional  analog  signal  processing,  GOMAC  03,  Tampa, 
FL,  Apr.  2003. 

13.  Kelvin  H.  Wagner,  Zoya  Popovic,  Dana  Z.  Anderson,  Randall  W,  Babbitt,  Lloyd  Grif¬ 
fiths,  Andre  Knoesen,  “RF  Photonic  Systems  for  Array  Control  and  Processing,”  In¬ 
vited  MURI  overview  at  GOMAC,  Monterey  CA,  March  2002 

14.  F.  Schlottau  and  K.  Wagner,  “RF  Photonics  for  Simultaneous  Multiple  TTD  Beam¬ 
forming  for  2-D  Antenna  Arrays,”  The  International  Topical  Meeting  on  Microwave 
Photonics  MWP  2001,  January  2002. 

15.  K.  Wagner,  Greg  Kriehn,  Friso  Schlottau,  Wideband  All-optical  BEAMTAP,  The  In¬ 
ternational  Topical  Meeting  on  Microwave  Photonics  MWP  2001,  January  2002. 

16.  G.  Kriehn,  F.  Schlottau,  and  K.  Wagner,  “Optically-Implemented  2-D  Beam  Steering 
and  Jammer  Nulling  using  BEAMTAP,”  submitted  to  ICO  Topical  Meeting  on  Optics 
in  Computing,  SRIE,  April  2002. 

17.  Kelvin  H.  Wagner,  Friso  Schlottau,  and  Jaap  Bregman,  Array  Imaging  Using  Spatial- 
Spectral  Holography,  submitted  to  ICO  Topical  Meeting  on  Optics  in  Computing, 
SRIE,  April  2002. 

18.  J.A.  Hagerty,  Z.  Popovic,  Passive  Millimeter- Wave  Ranging  using  Discrete  Lenses  with 
Wave- Front  Coding,  2001  European  Microwave  Conference  Digest,  pp.421-424,  Lon¬ 
don,  Sept.  2001. 

19.  J.  Vian,  Z.  Popovic,  “Smart  lens  antenna  arrays,”  2001  IEEE  International  Microwave 
Symposium  Digest,  pp.129-132,  Phoenix,  Arizona,  2001. 

20.  D,v  Anderson,  E.  Fotheringham,  S.  Romisch,  P.  Smith,  Z.  Popovic,  “Smart  antenna 
arrays  with  adaptive  optical  processing,”  IEEE  Trans.  Antennas  and  Propagation, 
Special  Issue  on  Wireless  Communications,  pp.  607-617,  May  2002. 

21.  P.  Smith,  E.  Fotheringham,  D.  Anderson,  Z.  Popovic,  “Adaptive  signal  processing  for 
microwave  carrier  broadband  signals,”  2003  GOMAC  conference  digest,  Tampa,  FL, 
April  2003. 

22.  R.  Reibel,  Z.  Barber,  M.  Tian,  W.  R.  Babbitt,  ’’Temporally  Overlapped  Linear  Fre¬ 
quency  Chirp  Programming  for  True  Time  Delay  Applications,”  poster,  Holeburning, 
Single  Molecules,  and  Related  Spectroscopies,  November  18-23,  2001,  Taipei,  Taiwan. 

23.  M.  Tian,  R.  Reibel,  Z.  Barber,  W.  R.  Babbitt,  “Broadband  true-time  delay  in  Tm:YAG,” 
poster,  Holeburing,  Single  Molecules,  and  Related  Spectroscopies,  November  18-23, 
2001,  Taipei,  Taiwan. 

24.  M.  Tian,  R.  Reibel,  Wm.  R.  Babbitt,  “Demonstration  of  broadband  true-time  delay 
with  optical  coherent  transient,”  MWP’2001,  January  2002,  Long  beach,  CA. 

25.  M.  Tian,  Z.  Barber,  T.  Chang,  R.  R.  Reibel,  and  W.  R.  Babbitt*  (*  Invited  Speaker), 
“The  Effects  of  Optical  Nutation  on  Stimulated  Photon  Echoes,”  Holeburning,  Single 
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Molecules,  and  Related  Spectroscopies,  November  18-23,  2001,  Taipei,  Taiwan.  (Re¬ 
lated  work  funded  by  AFOSR) 

26.  M.  Tian*,  R.  R.  Reibel,  Z.  Barber,  J.  Fischer,  W.  R.  Babbitt  (*  Invited  speaker), 

Optical  Coherent  Transient  True-time  delay:  Broadband  programming  methods,” 
2002  Physics  of  Quantum  Electronics  Conference,  Snowbird,  UT,  January  2002. 

27.  Joe  Fischer,  Zeb  Barber,  Randy,  Reibel,  Mingzhen  Tian,  and  Randy  Babbitt,  ” Linear 
Phase  Chirp  Programming  for  OCT’s,”  Optec2001,  (Bozeman,  MT,  August,  2001) 

28.  M.  Tian,  R.  R.  Reibel,  Z.  Barber,  J.  Fischer,  and  W.  R.  Babbitt,  “Optical  Coherent 
Transient  True-Time  Delay  Generator  Programmed  with  Linear  Phase  Chirp,”  in  Con¬ 
ference  on  Lasers  and  Electro-Optics,  Technical  Digest  (Optical  Society  of  America, 
Washington  DC,  2002),  May  2002,  Long  Beach,  CA. 

2.6.3  Patents  and  Applications 

1.  Invention  disclosure  submitted:  P.  Smith,  D.  Anderson,  Z.  Popovic,  E.  Fotheringham, 

A  dynamic  system  implementation  of  analog  optical  independent  component  analy¬ 
sis.” 

2.  K.  D.  Merkel,  Z.  Cole,  K.  Rupavaratham,  W.  R.  Babbit  and  K.  H.  Wagner,  US  Patent 
Application  No.  PCT/US03/14,612,  “Method  and  Apparatus  for  Processing  High 
Time-Bandwidth  Signals  Using  a  Material  with  Inhomogeneously  Broadened  Absorp¬ 
tion  Spectrum.” 


2.7  Classes  taught  by  Pis  relating  to  MURI  research 

1.  Dr.  Wagner  and  Dr.  Anderson  jointly  taught  a  graduate-level  Advanced  Optics  Lab 
in  Fall  of  2001  developed  for  interdisciplinary  students  with  interests  in  Chemistry, 
Physics,  and  Electrical  Engineering.  Several  of  the  MURI  graduate  students  were 
involved  either  as  TAs  or  in  taking  the  class. 

2.  Zoya  Popovic  taught  a  class  in  RF-optical  techniques  at  the  Technical  University  in 
Munich  in  2001.  This  was  a  part  of  an  ongoing  collaboration.  She  will  teach  again  in 
the  summer  of  2003. 

3.  Kelvin  Wagner  taught  CU  ECEN  6006,  Applied  Photonics,  new  course  developed  to 
be  the  second  course  in  the  new  core  optics  curriculum.  Spring  2002. 

4.  Kelvin  Wagner  taught  CU  ECEN  4016/5016,  Fundamentals  of  Photonics,  new  intro¬ 
ductory  course  developed  as  an  introduction  to  the  optics  core  curriculum.  This  was 
the  second  offering  of  this  course,  and  the  first  by  Prof.  Wagner. 
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5.  Randy  Babbitt  taught  MSU  Physics  353,  Holography-Laser  Photography,  Spring  2001, 
undergraduate.  Spring  01  (Randy  Babbitt),  Spring  02  (Randy  Reibel  taugth  it,  Randy 
Babbitt  was  the  mentor). 


Optoelectronic  Signal  Processing  Using  Finite 
Impulse  Response  Neural  Networks 
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work  in  the  above  mentioned  discipline. 


Xavier  da  Silveira,  Paulo  Eduardo  H.  B.  (Ph.D.,  Electrical  Engineering)  number  of  jammers  and  the  final  nulling  depth.  Two  sets  of  simulations  are  presented: 

Optoelectronic  Signal  Processing  Using  Finite  Impulse  Response  Neural  Networks  fc^e  set  applied  to  RF  antenna  arrays  and  the  second  set  applied  to  an  experimental 

Thesis  directed  by  Prof.  Kelvin  H.  Wagner  implementation  of  a  sonar  adaptive  array. 

Experimental  results  are  presented  for  a  single-layer  optical  processor  making 
This  thesis  investigates  the  use  of  finite  impulse  response  neural  network  as  the  use  of  a  scrolling  spatial  light  modulator  for  representing  the  input  signal  and  its  de- 
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application:  the  processing  of  signals  from  adaptive  antenna  arrays.  This  research  is 
initiated  by  computer  simulations  of  different  scenarios  with  multiple  broadband  sig¬ 
nals  and  jammers,  in  planar  and  circular  arrays,  studying  issues  such  as  the  effect  of 
modulator  non-linearities  to  the  performance  of  the  array,  and  the  relation  between  the 
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A  Photorefractive  grating  formation  and  diffraction  in  FIRNN  architectures 
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4.63  Receptivity  pattern  of  a  circular  array  in  the  presence  of  a  repeating 

broadband  jammer.  The  jammer  is  linearly  transformed  into  the  desired 
signal  instead  of  being  nulled . 

4.64  Learning  curve  in  the  presence  of  a  broadband  jammer  in  a  circular  array 

using  an  input  matrix  with  2048  time  samples . 


5.1  Optoelectronic  architecture  for  a  sonar  adaptive  array  processor  using  5.4  Schematic  diagram  of  the  digital  control  section  of  the  experimental 

a  scrolling  SLM  for  delayed  input  representation.  Dynamic  gratings  setup.  The  SLM  driver  provides  the  SYNC  signal  used  to  synchronize  the 

recorded  in  a  volume  hologram  are  used  for  implementing  the  adaptive  in-  HC11  micro-controller.  The  output  signal  0(k)  is  digitized  at  every  time 

terconnections,  followed  by  spatial  integration  onto  the  output  detectors.  step  during  the  read  cycle.  It  is  subtracted  from  the  desired  signal  stored 

Differential  heterodyne  detection  is  used  for  bipolar  signal  representa-  in  its  internal  EEPROM  and  the  resulting  error  signal  E(k)  is  stored  in 
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an  optic  beam  at  its  first  order  of  diffraction  modulated  by  the  error 
signal  e(«)  at  RF  carrier  frequency  a/i,  plus  a  RF  tone  at  o;2.  The  second 
AOM  diffracts  the  same  beam  at  its  -1  order,  producing  a  final  signal 
modulated  by  e(«)  at  RF  baseband,  plus  a  tone  at  the  IF  u\  -  w2-  .  .  . 
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5. 12  Anti-domain  visualization  in  a  SBN  crystal  before  the  repoling  attempt,  showing  the  detection  of  a  broadband  signal  in  the  frequency  range  of  the 

showing  a  high  concentration  of  180°  domains .  163  desired  chirp.  The  dashed  lines  depict  the  initial  and  final  frequencies  of 

5.13  Anti-domain  visualization  after  the  repoling  attempt,  showing  almost  no  the  original  desired  chirp . 

modification  in  the  180°  domain  distribution .  163  5.21  SLM  image  diffracted  from  the  PRC,  showing  the  correlation  pattern  of 


5.29  The  Fourier  plane  of  the  input  scrolling  SLM  for  different  jammer/desired  5.40  One  of  the  128  images  captured  from  the  SLM,  showing  a  broadband 

signal  power  ratios  (assuming  jammer  power  »  noise  power) . 180  chirp  and  a  broadband  jammer  detected  by  a  circular  array. . 

5.30  CW  jammer  nulling  and  near-field  beam-forming:  spectra  of  the  output  5.41  Image  of  the  SLM  in  the  Fourier  plane,  showing  a  uniform  signal  distri- 

of  the  processor  for  input  SINRs  ranging  from  0  dB  to  -6.02  dB . 181  bution  over  all  possible  transverse  spatial  frequencies,  characteristic  of  a 

circular  array. . 
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in  order  to  achieve  strong  non-linearities  without  becoming  a  bottleneck  to  the  flow  of 
information.  Such  an  architecture  is  commonly  found  in  neural  networks,  where  the 
nodal  interconnections  usually  require  a  high  order  of  linear  operations  (N2  for  inter¬ 
connecting  N  nodes  to  N  nodes)  followed  by  a  lower  order  of  non-linear  operations  (AT), 
as  depicted  in  figure  1.1. 
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matrix,  producing  all  inner  products  between  the  two  matrix  elements.  These  products  a  plane  WaVe  reference  beam  wbicb  18  focused  onto  the  CCD.  The  tilt  of  the  reference 

are  then  added  in  an  electrical/parallel  circuit,  producing  the  necessary  matrix  multi-  beam  can  be  varied  through  the  rotatinS  mirror  (RM),  such  that  for  each  input  pattern  a 

plication  between  the  input  and  weight  matrices.  Learning  is  performed  in  a  computer,  different  tilt  is  applied  to  the  mirror,  resulting  in  a  focal  spot  at  different  spatial  positions 

where  the  output  Y  is  compared  with  a  target  input  Y,  producing  an  error  term.  The  at  th®  CCD'  ThuS’  th®  grating  in  PRC2  maps  the  input  traininS  set  distinct  spatial 

spots  at  the  output  CCD. 


r 


a  g 

a  £ 

£  CO 

JS  1 

+3  £ 


S3 

P, 

443 

S3 

O 


2 

<*4-H 

I 

.© 

*-(-< 

S 


ss 

I 

£ 


•a  'S 


c$ 

43 


!§ 

d 

o 

o 

-u 

a 

4-3 

S3 

O 

,8 


o 

© 

I 

1 


£ 

l* 

& 

ctf 

3 

c 

03 

S3 

d 

| 

9- 

■s 

Cm 

o 

3 

0 


CD 


lO 

.CO. 

■g 

o 

a 

o 


♦s*  *-i 


■s 

& 

© 

o 


a 

o 

•8 


hC 

S3 


SP 


*  -3 


•S' 

1 


I 

43 

4J 

* 


i  1 


a 

J 

O 

*3 

8 

-+J 

O 

-2 

0) 

o 

-M 

a 

O 


CO 


a 

o 


& 


a 

a 


a 

4h 

o 

d 

43 

I 

o 

'o 

a 

4-3 

§ 

*— < 

g 

CJ3 

bp 

1 

43 

>> 

1 

.3 

a 

a 

*4 

g| 

s 

*3 

*9 

u 

O 

a 

•S 

4-3 

S3 

bfl 

•& 

jt 

a 

“© 

© 

© 

§ 

1 

4-3 

i 

a 

8 

bJO 

.3 

*53 

© 

§ 

i 

.s 

CO 

8 

a 

S3 

CO 

CO 

2 

> 

© 


I 


«3 

S3 

S3 

! 


xi 

© 


43 

I 

Js 


J 


o 

£ 


43 

H 

4*5 

t-H 

1 

© 

S3 


5 


CO 

co 


S’ 

*3 

c3 

© 


*3 

5 

2 

•a 


& 

Js 

© 

£ 

3 


i 


s? 


© 

a 

3 

© 

I 


S3 
o 

22  «s 

«2  ‘JH 

a  ° 

J?  © 

C  ^ 

X 

9 

§ 

o 


3° 

ud 


•8 

-w 

o 

a 

© 


1 


£ 


§ 

4-3 

*§ 

a 

bO 

.3 

a 

4-3 

o 


"S 


o 

bp 


P? 


fl3 

_© 


3 

*3 

o 


JL> 

3 


a 

4-3 

CO 

g* 

i 

4-3 

axis, 

<1 

£ 

-g 

a 

fe 

© 

a 

a 

% 

© 

43 

S3 

.2 

a 

© 

443 

443 

d 

a 

1 

1 

-0 

bP 

a 

O 

443 

443 

o 

© 

a 

S 

a 

g 

© 

-2 

a 

© 

43 

443 

*3 

8 

a 

(4 

o 

a 

a 

s 

a 

4-3 

.S5 

43 

43 

4-3 

g 

d 

.3 

1 

4-3 

43 

S3 

O 

d 

u 

fcuO 

*4 

d 

i 

t6 

4-3 

44) 

a 

I 

i 

d 

CO 

© 

1 

© 

.3 

8 

& 

.& 

a 

4-3 

| 

443 

.2 


§ 


1 


s  g 

S  “ 


3 

S3 


H 


d 

t: 

I 

s 

t9 

J§ 


o 

S3 

-Q 


§ 


a 

H 

S3 

.2 

*«5 

S3 

a 


S3 

O 


S3 

O 


I  1 


& 

© 


ti  .3 


§ 


•§ 


S3 

a 

.3 

43 

4-3 

43 

4-3 


fab 

■2 


! 

CO 

S 

1 

bC 

a 

J 

© 

43 


O 

a 

o 


3  5 

a  *s 


o 

*> 

© 


© 


Q 

a 


d 

'5s 

o 

© 


45 

a 


*5* 

8 

.3 

a 

4-= 

& 

© 

ti 

a 


© 

45 


S? 


hO 

© 


hO 

.3 

"3 

fcb 

a 


& 

9 


43 

H 

d 

o 

o 

4-3 

S3 

a 

4-3 

S3 

O 

J§ 


o 

a 

a 


o 

a 


© 

a 


d 

a 

4-3 

S3 

a 

S3 


■8 


3 


I 


,2  43 


a 

a 

1 


© 

43 


d 

a 


S3 

8 

,© 

53 

x 

d 

O 

a 

o 

■8 

bjO 


C? 

"a 

a 


d 

a 


ceivers,  where  an  adaptive  algorithm  steers  the  antenna  gain  pattern  towards  sources 
that  are  coherent  with  the  desired  signal.  The  same  antenna  arrays  that  are  used  for 
reception  can  also  be  used  for  transmission.  If  the  weights  are  kept  constant,  the  di¬ 
rectivity  pattern  of  the  array  used  as  a  transmitter  will  be  the  same  as  the  receptivity 
pattern  of  the  same  array  used  as  a  receiver. 


t3  .W 


£  TJ 


T  taps  for  every  antenna  element.  Such  a  TDL  array  can  form  N  broadband  beams 
nulls  simultaneously  without  beam  squint. 
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the  previous  layers  are  multiplied  by  the  weights,  summed  and  delayed,  evaluating  temporal  correlation.  The  output  of  each  FIR  is  then  multiplied  by  a  delayed  version  of 

equation  2.2.  Simultaneously,  a  derivative  of  the  non-linear  function  with  respect  to  non"^near  derivative,  producing  the  delta  terms  used  to  update  the  previous  layer, 

the  partial  sums  is  generated.  The  output  error  is  calculated  at  the  output  of  the  we*Shts  are  then  updated  by  an  accumulation  of  the  outer  product  between  the 

network.  This  term  is  multiplied  by  the  learning  rate  and  the  non-linear  derivative,  back-propagated  S  terms  and  a  delayed  version  of  the  forward  propagating  input  signal, 

producing  the  delta  terms  used  to  update  the  weights  in  the  output  layer.  Part  (b)  There  are  three  important  points  to  be  observed: 


temporal  dot-product  found  in  equation  2.22  and,  therefore,  the  same  time  and  space- 
integrating  methods  can  be  used,  followed  by  an  inner-product  with  vector 

rn)). 


m 

in 


ing  operations  required  in  the  implementation  of  FIRNNs,  and  their  efficient  imple¬ 
mentation  is  critical  in  forward  and  back-propagation,  in  the  conventional  and  in  the 
delayed-feedback  algorithms. 


containing  the  signal  £(«)  is  set  to  interfere  with  the  reference  beam.  The  two  beams  limitation  can  be  compensated  by  using  a  high  resolution  device.  Figure  3.3  shows  the 

continually  record  holographic  gratings  in  the  crystal  as  it  rotates  at  an  angular  velocity  resulting  image  when  a  256  x  256  FLC-SLM  is  fed  with  a  binary  chirp  with  frequency 

0.  Since  the  reference  beam  is  collinear  with  the  axis  of  rotation,  it  is  always  Bragg  varying  linearly  from  93.8  Hz  to  563  Hz  at  a  1.50  kHz  sampling  frequency.  All  input 

matched  with  the  previously  recorded  gratings.  As  a  consequence,  the  delayed  signal  nodes  receive  the  same  signal,  but  with  a  delay  of  0.18  ms  between  nodes,  which  causes 

0  is  spatially  mapped  into  a  cone  of  simultaneously  diffracted  holograms,  where  a  sPa^ai  shift  of  the  resulting  pattern.  Figure  3.4  shows  the  Fourier  plane  of  this 


Holographic  weights 


grating,  decreasing  its  maximum  dynamic  range. 


generate  the  weight  update  term  £(2)(«)  that  is  fed  to  the  1-D  SLM  4.  This  plane  pixel  aligned  through  PRC  1  onto  detector  array  1.  The  first  layer  requires  the  use  of  a 

is  imaged  through  PRC  2  onto  detector  array  2.  The  second  scrolling  SLM  displays  delayed  version  of  its  input  for  updating  the  weights,  since  equation  2.16  in  the  2-layer 

2<1>(k).  Interference  of  this  signal  array  with  the  error  gradient  signal  coming  from  case  becomes  A w$r  =  -»tfjl)(#c  -  T)s[0){k  -  T  -  r).  This  can  be  achieved  in  two 

SLM  4  produces  the  outer-product  used  for  updating  the  weights  in  the  volume  hologram  different  ways:  the  first  one  consists  of  using  a  periodic  input  signal  with  a  period  of  T 

PRC  2.  The  grating  evolution  in  PRC  2  is  given  by  ^  time-steps,  and  using  an  input  SLM  with  T  rows.  In  this  way,  s|0)(k-t)  =  s^(k-T-t) 
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integration  of  the  partial  products,  while  the  second  performs  the  interconnections  first  create  hybrid  architectures  by  freely  combining  layers  implemented  by  the  conventional 

followed  by  spatial  integration  and  time  integration  in  a  TDI  photodetector.  temporal  back-propagation  architecture  with  other  layers  implemented  by  the  delayed 

These  two  basic  architectures  were  extended  and  provided  with  closed-loop  feed-  feedback  architecture.  This  is  possible  since  the  electronic  signals  3*(/c)  and  S1+\k  - 

back  in  order  to  implement  FIRNN  architectures  capable  of  on-line  gradient  descent.  T^),  the  only  ones  passed  between  layers  l  and  l  +  1,  are  functionally  identical,  no 

Two  single-layer  architectures  were  presented:  one  based  on  the  input  delay  plane  ar-  matter  which  architecture  is  being  used. 


while  the  linear  order  non-linear  operations  have  been  performed  electronically,  thus 
exploiting  the  potential  advantages  of  each  technology. 

Finally,  the  2-layer  architectures  presented  could  potentially  be  extended  to  an 
arbitrary  number  of  layers,  where  special  care  should  be  taken  in  the  amount  of  delay 
applied  to  the  fixed  delay  lines.  Note  that  in  a  multi-layer  implementation  one  could 
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The  software  used  in  the  generation  of  all  of  the  simulations  presented  in  this 
chapter  is  explained  in  detail  in  appendix  B.  This  chapter  is  divided  in  two  main 
sections:  Section  4.2  describes  simulations  conducted  using  the  BEAMTAP  version  of 


antenna  array  processor  requiring  the  use  of  only  two  tapped  delay-lines.  Figure  4.1  partial  products  are  temporally  delayed  and  accumulated  within  the  output  delay-line 

illustrates  the  simulation  of  this  architecture  by  showing  some  of  the  signals  captured  (top  of  FiS*  4*1)»  producing  an  output  signal  o(«)  which  is  subtracted  from  the  desired 

during  a  snapshot  of  the  simulation  during  the  processing  of  the  signals  in  a  64-element  signal  d(«),  producing  the  feedback  error  signal  e(/c).  This  signal  is  fed  into  the  tap-out 

array  with  64  taps  in  its  two  delay-lines.  The  leftmost  graphic  depicts  the  spatio-  delay-line  and  its  outer  product  with  a  delayed  version  of  the  input  [s(k  —  T)]  is  used  to 

temporal  RF  field  amplitude  sampled  at  8  GSamples/sec  along  the  horizontal  time  adapt  the  weights.  At  steady  state,  the  weights  accumulate  the  cross-correlation  of  the 


Desired  signal  ^formation  8^nat , 9Pectrum  Info  signal  and  error  corr. 


signal.  Figure  4.5  depicts  the  spectrum  of  an  unknown  PNBPSK  (pseudo- noise  bipolar  figure  4.7  shows  a  plot  of  the  cross-correlation  between  the  feedback  error  signal  and 

phase  shift  keying)  signal,  which  is  first  generated  by  a  pseudo-random  number  sequence  the  pseudo-random  sequence.  The  strong  correlation  peak  proves  that  the  pseudo- 

making  it  spread-spectrum),  and  then  modulated  to  the  same  center  frequency  and  random  sequence  is  present  in  the  error  signal  and  allows  for  various  spread-spectrum 

with  the  same  bandwidth  as  the  desired  signal  in  order  to  produce  an  overlap  of  their  communication  techniques  to  be  used  in  order  to  convey  unknown  information  to  the 


Frequency  response  at  different  angle*  Frequency  response  at  different  anglee 
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exactly  at  the  mid-point  of  the  sinusoidal,  even  orders  of  non-linearity  also  occur.  linear  terms,  as  will  be  shown  in  section  4.2.6.  Note  that  usually  only  the  positive 

When  a  temporal  signal  with  frequencies  components  f\  and  f<i  undergoes  non-  temporal  frequencies  are  shown  in  the  simulations,  but  their  negative  conjugates  are 

linear  modulation,  the  even-ordered  non-linear  terms  give  rise  to  an  output  signal  of  always  present  since  all  simulated  signals  are  real. 

frequency  components  2/i,  2/2,  4/i,  4/2,  f\  —  fi,  fi  —  /1,  f\  +  /2,  ...  The  odd-ordered  Figures  4.22  and  4.23  illustrate  this  principle  through  a  simulation.  The  first 

terms,  on  the  other  hand,  give  rise  to  the  frequency  components  3/i,  3 /2>  5/i,  5 /2,  figure  (Fig.  4.22)  shows  a  broadband  chirp  detected  in  the  far-field  with  an  AOA  of 


spatio-temporal  spectrum  of  the  desired  signal,  acting  as  an  in-band  main-beam  jam- 
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Figure  4.56  shows  the  corresponding  spatio-temporal  Fourier  spectrum  of  the  detected 


Learning  curve. 


(4)  A  linear  adaptive  array  is  capable  of  adapting  for  the  detrimental  effects  of 
non-linearities.  However,  the  available  number  of  degrees  of  freedom  is  reduced, 
and  one  can  devise  instances  in  which  the  performance  of  the  array  degrades 

rapidly  (such  as  main  beam  jamming).  Chapter 
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true-time-delay).  This  task  becomes  specially  difficult  when  the  signals  are  present  in 
the  near-field  of  the  array,  or  when  the  array  shape  is  not  planar  or  is  time  varying, 
as  these  scenarios  require  a  higher  number  of  degrees  of  freedom  and,  therefore,  more 
processing  capacity.  Several  photonic  implementation  schemes  have  been  previously 
proposed  [68,  70,  58]  and  discussed  in  chapter  1.  However,  we  believe  that  this  is  the 
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rising  edge  is  detected,  the  output  signal  0{k)  is  read  out  through  the  A/D  (analog-to-  software  package  [101].  Note,  however,  that  the  DAQ  board  and  the  PC  are  present 

digital)  input  port  E  (pin  1)  and  it  is  subtracted  from  the  desired  signal,  which  has  been  in  order  to  program  the  micro-controller  and  to  provide  a  second  way  for  collecting 

previously  programmed  in  the  internal  EEPROM  of  the  micro-controller.  The  resulting  output  data  (besides  probing  signal  0(k)  with  an  oscilloscope),  and  their  presence  is 

128  error  signals  E(k)  are  stored  in  RAM  until  the  next  write  cycle.  Its  mean  value  not  required  for  the  normal  operation  of  the  adaptive  processor, 

is  calculated  and  subtracted  from  each  sample  in  the  time-series,  in  order  to  compen-  The  D/A  (digital-to-analog)  converter  (DAC0800)  and  the  LM741  op-amp  con- 


150 

vert  the  digital  output  into  a  zero-balanced  analog  waveform.  The  LM6221  unit  gain 
buffers  [102]  drive  the  50  «  lines  used  to  drive  RF  mixer  before  the  RF  amplifiers  and 


two  stages,  but  it  is  high  enough  to  saturate  the  last  two.  Once  saturated,  op-amps  *  Pr°Vide  maXimUm  “"P15^6  at  the  DSB  modulated  si^  at  which 

require  some  time  to  return  to  normal  operation,  and  such  a  delay  translates  into  a  P°int  **  ^  difference  +•  ~  *  to  zero-  Passive  stabilization  guarantees  that 

degradation  of  the  temporally  multiplexed  signal  at  the  beginning  of  every  read  cycle.  thiS  Wi“  be  ^  for  *  feW  SeC°ndS  and  pr°videS  enough  time  t0  initiate  Mtive 

In  order  to  prevent  this  degradation,  the  CLC532  [107)  high  speed  analog  multiplexer  stabilization.  Active  heterodyne  stabilization  is  performed  using  a  SR850  DSP  Lock-in 

amplifier  in  combination  with  a  mirror  attached  to  a  piezo-electric  actuator,  as  shown  in 
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Figure  5.14  shows  a  map  of  the  birefringence  of  the  BaTiOj  crystal.  It  is  done  by 
Ail  the  following  experimental  results  reported  on  this  thesis  have  been  obtained  imaging  the  crystal  while  illuminated  by  a  single  circularly  polarized  beam  propagating 

using  a  45°-cut  ferroelectric  Fe-doped  barium  titanate  (BaTi03)  crystal  of  dimensions  normally  to  the  crystal  surface,  followed  by  an  analyzer  placed  after  the  crystal.  With 

[5x5x3]mm3,  which  also  has  a  4mm  point  group  symmetry  at  room  temperature.  The  a  perfectly  collimated  beam,  regions  with  the  same  amount  of  birefringence  change  the 

polarization  of  the  incident  beam  by  the  same  amount  and,  therefore,  present  equal 
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^  **4.  f  •  i.  OT  Ayr  r  u  .  -  .  ,  .  the  inPut  SINR  decreases,  evidenced  by  the  dimming  of  the  broadband  section  in  the 

,  producing  the  set  of  input  SLM  frames  shown  m  the  Fourier  plane  in 

.  .  .  .  ,  ,  Fourier  plane.  How  does  the  decrease  in  the  input  SINR  affect  the  performance  of  the 

that  the  jammer  becomes  increasingly  predominant  over  the  signal  as 

processor?  This  questions  is  answered  by  the  plots  in  figures  5.30  and  5.31,  which  show 


(gp)  j««od  jndino  jiaonsuj 


space,  the  desired  chirps  are  clearly  recovered  and  coherently  added  in  order  to  provide 

„  .  i  ,  .  A  ,  o  AOA  of  the  far-field  source  (see  Section  5.5  for  more  details).  Compare  this  fitrure  to 

a  strong  and  clean  output  signal.  The  BER  for  this  scenario  is  6.20xl0~3.  H  ngure  to 

n*  ,QQ.  ,,  .  , .  A  figure  4.57  and  note  the  remarkable  similarity  with  the  weights  after  adaptation  in  the 

Figure  5.38  shows  the  correlation  pattern  after  adaptation.  Note  that  two  sinu-  6  H 

. ,  ,  ...  ,  simulation  of  multiple  beam- forming  in  a  circular  array, 

soidal  correlation  patterns  are  present,  each  corresponding  to  a  response  towards  each 

•  j*  •  i  i  i*  mt  i  -  ,  . .  Figure  5.39  shows  the  correlation  pattern  in  the  Fourier  plane.  Compare  it  to 

one  of  the  individual  chirps.  The  phase  of  each  sinusoid  is  shifted  according  to  the  H  P 
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a  linear  transformation  which  will  make  the  jammer  resemble  the  desired  signal  as  Figure  5.44  shows  the  correlation  pattern  in  the  Fourier  plane.  Note  that,  as  in  the  pre¬ 
close  as  possible,  artificially  resulting  in  a  better  SINR  than  expected.  Such  an  effect  vlous  8cenario> the  processor  responds  to  a  broadband  signal  coming  from  all  transverse 

was  simulated  in  section  4.3.5  in  chapter  4,  with  similar  results.  The  BER  achieved  spatial  frequencies, 

in  this  scenario  is  6.02  xlO-3.  It  is  due  to  a  limitation  in  the  SLM,  only  capable  of 
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Figure  5.48  shows  the  correlation  pattern  after  adaptation.  Note  that  again  a 

thin,  curved  line  is  obtained,  as  expected  from  a  broadband  signal  in  the  near-field.  5‘4'8  Far'fleld  be*™-f°rming  and  multiple  jammer  nulling 

However,  the  shape  of  the  line  is  different  from  that  presented  in  figure  5.27,  since  In  the  finaJ  scenario,  the  processor  performs  beam-forming  in  the  direction  of 

the  source  is  present  at  a  different  range  and  AOA  this  time.  Figure  5.49  shows  the  a  broadband  chirp  in  the  far-field  of  the  antenna  array  (AOA  =  0.25  radians)  while 
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As  shown  in  the  previous  scenarios,  it  is  possible  to  obtain  a  pattern  proportional  ^  °  ^  the  intel"element  antenna  «  *  the  velocity  of  the  acoustic  wave  in 

to  the  correlation  between  the  desired  signal  and  the  error  signal.  To  a  first  order  ^  °  ^  **“  A°A  °f  the  >ignaI*  ThuS’  the  P°siti°n  °f  the  correlation  Peak 


(3)  Time-multiplexing  can  be  successfully  employed  for  holographic  read/ write  i: 
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one  based  on  the  delayed  feedback  algorithm,  known  as  the  Broadband  and  Efficient  presented  for  planar  and  circular  adaptive  arrays,  with  signals  present  in  the  far- field 

Adaptive  Method  for  True-time-delay  Array  Processing  (BEAMTAP)  optoelectronic  and  in  the  near-field.  Succesful  beam-forming  was  observed  in  all  scenarios.  Jammer 

architecture,  and  the  other  based  on  the  conventional  temporal  back-propagation  al-  nulling,  however,  was  only  effectively  observed  in  the  case  of  narrowband  jammers, 

gorithm,  known  as  the  Finite  Impulse  Response  Optical  Neural  Network  (FIRONN)  The  limited  memory  of  the  scrolling  SLM  driver  made  it  necessary  to  repeat  the  input 

architecture.  The  simulations  focused  on  the  processing  of  signals  from  adaptive  an-  frames  after  128  samples,  resulting  in  a  partial  correlation  between  the  desired  signal 


and  broadband  jammers,  allowing  the  processor  to  linearly  transform  the  jammers  into  (1)  Improved  SLM  driver:  A  simple  way  to  improve  the  performance  of  the 

the  desired  signal  rather  than  nulling  them.  experimental  demonstration  would  be  to  replace  the  current  SLM  driver  with 

An  experiment  was  performed  in  order  to  determine  the  relationship  between  the  an  improved  model  with  twice  as  much  memory,  allowing  the  use  of  256  taps 

intensity  of  a  CW  jammer  and  the  processing  gain  of  the  processor,  showing  that  the  instead  of  128.  Better  yet,  a  driver  that  could  accept  data  in  real-time  could 

processing  gain  increases  roughly  linearly  with  the  intensity  of  the  jammer.  Finally,  it  be  used.  This  would  overcome  the  problem  arising  from  the  partial  correlation 
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Let  and  EB(xo,yo;  t)  represent  the  time-varying  forward  propagating 

signals  whose  interference  pattern  records  a  dynamic  grating  in  the  PR  crystal.  Their 
modulating  SLMs  (SLM  1  and  SLM  2)  and  detectors  (Detector  array  1  and  Detector 
array  2)  are  one-dimensional  and  two-dimensional,  respectively.  Note  that  such  a  setup 
is  sufficient  to  describe  all  the  FIRNN  architectures  presented  in  chapter  3  (input  delay 


(3)  a  source  of  input  data  which  can  be  used  in  experimental  implementations,  such 
as  the  FIRONN  architecture  described  in  chapter  5. 


diagrams  illustrating  the  BEAMTAP  architecture. 
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from  the  main,  menu,  as  explained  in  section  B.2.  It  is  recommended  that  the  Array  and  •  X_TALK  —  the  amount  of  cross-talk  between  a  given  antenna  element  and  all 

Optic  parameters  be  defined  first,  followed  by  the  Jammer  and  Desired  signal  parameters  of  its  neighbors.  The  cross-talk  exponentially  decreases  with  physical  distance, 

and  the  Algorithmic  parameters,  in  that  order.  Each  one  of  these  set  of  parameters  is  Set  it  to  zero  for  no  cross-talk,  or  to  one  for  complete  cross-talk  (all  elements 

explained  next.  detect  the  same  signal). 


Spatial  appedization  factor  -  the  exponential  factor  used  in  the  spatial  •  Center  freq.  -  the  center  frequency  of  the  desired  signal. 

Gaussian  appedization  of  the  entire  array.  Set  it  to  zero  for  no  appedization. 

•  F  -  fractional  bandwidth:  the  ratio  between  the  total  signal  bandwidth  and  the 
rad  —  the  radius  of  the  circular  array.  Set  it  to  zero  to  define  a  planar  array.  desired  signal  center  frequency. 


•  Signal  wavelength  -  the  wavelength  of  the  center  frequency  of  the  desired  .  Vpi  -  voltage  V,  of  the  EO  modulator  (See  section  4.2.6  for  more  information), 

signal. 

•  bp  -  amplitude  of  uniformly  distributed  bias  voltage  variation  in  EO  modulators 
.  Element  spacing  -  the  spacing  between  array  elements,  set  in  order  to  Nyquist  across  antenna  elements.  Set  it  to  zero  to  prevent  bias  voltage  from  varying  from 

sample  the  input  signal.  it8  optimal  point> 


parameters  are:  to  accommodate  the  required  number  of  chirps  in  a  total  of  MAXJSTEPS 


gets  for  modifying  the  array,  optical,  input  signals  and  algorithmic  parameters, 
respectively. 
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This  thesis  entitled: 

Developments  in  Photorefractive  Two-Beam  Coupling  Systems 
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Department  of  Electrical  &  Computer  Engineering  The  final  copy  of  this  thesis  has  been  examined  by  the  signatories,  and  we 

find  that  both  the  content  and  the  form  meet  acceptable  presentation 

2000  standards  fo  scholarly  work  in  the  above  mentioned  discipline. 


and  oi  the  photorefractive  medium  in  a  single  element.  While  investigating  these 
spherical  crystals,  we  observed  internal  whispering-gallery  mode  oscillations.  These 
oscillating  patterns  adopted  simple  polygonal  shapes  depending  on  the  pump  entry 
and  crystal-axis  angles.  We  modeled  the  occurrence  conditions  of  the  lowest  order 


Interference  intensity  pattern 


polarization  state  of  the  fields.  In  this  treatment  we  only  take  into  account  the  main 
component  of  barium  titanate’s  electrooptic  tensor  (r42=1640  at  low  frequencies), 
since  the  other  components  are  more  than  one  order  of  magnitude  smaller.  For 


I 

1 

a 

a 

a 


for  the  signals  in  the  photorefractive  processors  is  given  by  Bragg  matching 
considerations  and  is  a  few  gigahertz  for  a  1  cm  grating.  This  limit  is  inversely 
proportional  to  the  grating  length. 


*3 


E 


entails  alignment  of  beams  through  crystals  or  other  optical  elements  which  are  part 
of  each  photorefractive  function.  This  means  that,  with  a  collection  of  appropriate 
modules  in  hand,  a  complex  system  can  be  quickly  assembled  or  reconfigured. 


terminated  optical  devices.  What  are  the  tradeoffs  of  such  benefits?  This  brief  section 
clarifies  the  way  in  which  the  use  of  multimode  fibers  affects  the  modules. 
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attached  with  epoxy,  and  its  end  polished  flush  with  the  ferrule's  face  to  provide  an 
optical  quality  surface.  The  prepared  ferrule  is  then  inserted  into  the  collet  and  a 
specially  designed  tool,  which  moves  the  collet-lens  assembly  with  respect  to  the 


didn’t  break.  And  most  importantly,  such  stress  is  very  likely  the  culprit  for  the  adjustment  range.  The  ellipses  represent  milled  slots  resulting  on  wide 

but  thin  (about  750  urn)  flexures. 

alignment  drift  as  the  aluminum  tries  to  relax. 


Figure  2.5.  Details  of  the  base  design.  Dimensions  are  in  millimeters, 
(a)  Original  design  with  vertical  pivots,  (b)  New  base  design  is  similar 
to  the  original,  but  with  horizontal  pivots  which  allows  for  greater 
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inside  the  crystal,  resulting  in  beam  diameters  in  the  waist  of  about  100  pm.  The  time  which  helps  isolate  them  from  external  vibrations.  The  buffered  fiber  also  has  the 

response  of  a  crystal  is  roughly  inversely  proportional  to  the  total  intensity  of  the  advantage  of  greater  mechanical  resistance,  are  white,  and  therefore  are  easy  to  see  in 

beams.  [Horowitz, '9 lb]  This  means  that  the  smaller  beam  diameters  should  result  in  a  a  darkened  room.  Not  only  did  we  use  these  buffered  fibers  in  some  of  our  new 

decrease  in  the  response  time  of  the  modules. 


AN  APPLICATION  EXAMPLE:  THE  AUTO-TUNING  FILTER 


experiences  greater  round-trip  loss,  which  results  in  an  even  weaker  signal,  and  so  on, 
while  the  stronger  signal  undergoes  the  opposite  behavior.  Ideally,  this  positive 
feedback  process  continues  until  the  weaker  signal  is  totally  inhibited  by  the  stronger 
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crystal  axis  continuously  and  permitted  the  angle  between  the  beams  to  be  chosen  3.4.  WHISPERING-GALLERY  MODE  PATTERNS 

from  5°  to  45°  in  5°  increments.  Maximum  gain  occurred  at  the  expected  45° 


Figure  3.2.  Two-beam  coupling  between  two  pairs  of  fibers  using 
spherical  photorefractive  crystal  of  barium  titanate. 


vertically  downward  direction.  The  pump  beam  is  focused  by  a  microscope  objective 
so  that  a  beam  waist  occurs  inside  the  disk  and  is  roughly  matched  in  position  with  a 
waist  of  the  mode.  Most  of  the  pump  energy  is  depleted  in  the  first  millimeter  after 
the  pump  enters  the  disk,  indicating  efficient  energy  transfer.  Back-side  oscillation 


As  further  investigation  we  calculate  the  gain  maxima  for  triangle  excitation 


also  observed. 


1  Phase-conjugation  observations  were  done  by  Dr.  Dana  Anderson  and  Dr.  Vladimir  Shkunov. 


opposed  to  the  techniques  mentioned  above,  is  adaptive  due  to  the  real-time 


optic  links  is  that  of  using  stimulated  Brillouin  scattering,  where  the  carrier  Different  versions  of  the  auto-tuning  filter,  some  with  fiber  optics  and  spherical  disks 

propagating  inside  an  optical  fiber  pumps  a  beam  which  propagates  in  the  opposite  like  the  ones  described  in  the  previous  chapter,  and  others  with  Gaussian  beam  optics 

direction.  Optical  feedback  of  the  generated  beam  enhances  the  carrier  suppression  and  other  crystal  cuts  have  been  tested.  However,  their  use  doesn’t  pertain  to  the 

and,  in  fact,  a  record  carrier  suppression  of  55  dB  has  been  reported  by  using  a  carrier  suppressor  and  will  not  be  considered  further  here. 


application.  We  first  discuss  how  the  performance  of  the  carrier  suppressor  described 
in  section  4.3  fulfills  the  performance  specifications  for  that  system.  We  then 
introduce  a  new  carrier  suppressor  which  takes  advantage  of  a  channelized 
electrooptic  modulator  design,  discuss  it’s  implementation,  and  expose  an  important 
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This  particular  set  of  data  was  taken  by  Amy  Sullivan. 


From  the  curve  in  Figure  4.8  we  obtain  the  curve  for  the  output  vs.  input 
modulation  strength,  which  is  shown  in  figure  4.9  (compare  to  the  theoretical  plots 

shown  in  figure  4.4).  This  curve  results  in  a  peak  output  modulation  strength  of  about  4.4.  FUTURE  WORK 

m(zh 44  dB. 


requirements,  and  discuss  its  current  limitations. 


Adequate  carrier  suppression  performance  is  demonstrated  in  this  chapter, 
however  the  incorporation  of  the  carrier  suppressor  in  the  antenna  array  system 
requires  a  design  that  can  process  several  channels  in  a  compact  space.  The  antenna 


On 


Dr.  Stefania  Romisch  collaborated  with  the  modulator’s  assembly  and  characterization. 


cause  is  the  RF  coupling  between  the  resonant  circuits  that  drive  the  electrooptic 


differential  equation  would  result  for  each  pair  of  plane  waves,  resulting  in  a 
complicated  problem.  With  the  operator  theory  the  spatial  basis  may  be  chosen  to  be 
a  set  of  orthogonal  facial  features,  bridging  the  gap  between  the  mathematical 
notation  and  the  information  processing  task  at  hand. 


[5.2].  The  field  values  at  any  position  z  within  the  crystal  is  determined  by  solving  _ 

J  y  ng  which,  by  eliminating  Ew(0),  gives  us  a  differential  equation  for  the  coupling 

for  the  coupling  operator  T.  As  we  will  see  later,  T  corresponds  to  a  rotation  operator 

operator: 

which  conserves  the  length  of  the  fieldvectors  (meaning  that  the  transformation  is 


For  this  simple  case,  closed-form  solutions  exist  for  both  the  density  operator  p  and  permutation  of  ij,k  and  -1  for  an  odd  permutation.  Third,  they  anti-commute, 

the  coupling  operator  T.  From  this  point  forward,  we  will  cease  to  refer  to  them  as  meaning  that  {o,.,o  j}  =  alaj  +<r p,  =  0  where  i*j. 

operators  but  instead  refer  to  their  2-by-2  matrix  representation.  The  coefficients  of  the  spin  matrices  can  be  conveniently  expressed  in 

Any  2-by-2  matrix  can  be  represented  by  a  superposition  of  the  unit  matrix  spherical  coordinates: 


As  can  be  deduced  from  the  expressions  above,  are  orthogonal  and  are  obtained  by 

rotating  c,  and  c2  around  the  o3  “axis”,  such  that  <s„  is  aligned  with  the  projection  of  =a3^20Oi  =G „  sin(29)  +o  3cos(20), 


5.2  (b)  shows  the  north  pole,  which  points  in  the  direction  of  a}.  Given  the  initial 


(a),  plotted  for  y=0.  Conversely,  an  imaginary  coupling  constant,  which  occurs  for 
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figure  5.2.  We  know  from  photorefractive  physics  that  a  real  coupling  constant, 
which  occurs  for  y=0  (see  equation  [5.2]),  corresponds  to  pure  energy  coupling  and 


The  fiber-interconnected  two-beam  coupling  modules  are  pre-aligned  building 
blocks  for  photorefractive  system  design.  They  provide  a  way  to  quickly  test  an  idea 
without  having  to  invest  the  time  in  critical  optical  alignment  which  is  necessary 
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Phys.  Rev.  A,  1995. 52:  p.  878-81. 

[Zozulya, '95].  Zozulya,  A.A.,  Saffman,  M.,  and  Anderson,  D.Z.,  Stability 
analysis  of  two  photorefractive  ring  resonator  circuits:  the  flip-flop  and  the 
feature  extractor.  J.  Opt.  Soc.  Am.  B,  1995. 12:  p.  1036-47. 
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An  Optically-Switched  Tlransmit/Receive  Lens  Array  for  Beam-Space  Adapti 
Communication  Systems 
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optically  controlled  transmit/receive  lens  array  is  developed.  The  optical  control  of 
the  lens  array  allows  the  array  to  switch  between  transmit  and  receive  modes  rapidly 
with  negligible  interference  to  the  microwave  signals.  A  low  optical  power  single  pole 
double  throw  switch  is  developed  for  routing  the  transmit  and  receive  signals  in  the 
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the  optical  fibers .  50  “b”  (red)  and  circuit  “e”  (green) . 

2.36  Illustration  of  the  test  setup  used  to  measure  the  unit  cell.  The  3.8  Measured  loss  in  circuit  “b”  for  the  switch  “on”  (red)  and  the 

test  setup  is  calibrated  to  an  aperture  the  size  of  the  unit  cell.  ...  50  switch  “off”  (blue) . 


4.4  Photograph  showing  an  aluminum  plate  that  allows  wire  bonding  sion  for  the  0°  off  axis  (a)  and  the  30°  off  axis  (b)  feeds.  The 

on  both  side  of  the  lens  array. . .  85  transmit  patterns  are  in  (blue)  and  the  receive  patterns  are  in  (red).  94 

4.5  Photograph  showing  the  active  side  (a)  and  passive  side  (b)  of  the 

RF  mount  for  the  unit  cell  with  their  alignment  pins .  86 


numbered  data  points .  174 

7.11  The  radiation  patterns  for  lens  array  system  in  multi-user  environ¬ 
ment  with  “desired”  user  at  30°  and  the  “other”  user  at  -40°  (a),0° 

(b),  20°  (c)  and  25°  (d) .  175 


The  concept  of  diversity  dates  back  to  the  1940's  when  it  was  used  to 
improve  radar  performance.  Much  of  the  early  published  work  focuses  on  diversity 
CHAPTER  1  reception  where  multi-path  fading  is  a  problem  [8,  9,  10,  11,  12].  For  example, 

transmission  over  long  distances  (i.e.,  part- way  around  the  world)  results  in  multiple 
TROnurTTOM  reflections  of  the  waves  off  the  ionosphere,  and  the  different  reflected  waves  interfere 


The  magnitude  and  phase  changes  in  the  received  signals  are  often  referred  significant  computational  power  to  calculate  the  adaptive  weights.  The  number 

to  as  weights,  in  reference  to  a  complex  weighted  sum.  The  weights  in  the  algorithm  of  computations  increases  geometrically  with  the  number  of  array  elements.  In  this 

have  an  interpretation  in  microwave  antenna  array  theory.  The  weights  are  the  thesis,  a  new  type  of  front  end  (antenna  array  and  receiver)  is  explored  in  the  context 

relative  magnitudes  and  phases  of  the  generators  at  each  antenna  element  for  a  of  adaptive  systems.  In  a  “standard”  planar  array,  each  element  receives  a  phase- 

transmit  array.  Because  of  reciprocity,  this  is  equivalent  to  the  amplitude  and  phases  shifted  version  of  the  same  input  plane  wave.  In  the  lens  array  presented  in  this  thesis, 


the  antenna  array  performs  a  Fourier  transform  and  images  the  signal  sources  onto  a 
focal  surface,  and  the  signal  image  is  then  sampled  with  receivers.  Most  applications 
require  both  transmit  and  receive  functions,  and  a  part  of  this  thesis  is  concerned 
with  fast  efficient  routing  of  signals  through  transmit  and  receive  paths.  To  that  end, 
the  transmit  and  receive  array  presented  in  this  thesis  uses  optical  control  which  does 
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The  adaptive  system  in  Figure  7.3  is  intended  for  a  narrowband  application 
of  bandwidth  BW,  where  the  delay  of  the  received  signals  across  the  array  or  the 
difference  in  delay  between  multi-path  signals  is  less  than  For  very  broadband 
systems,  such  as  some  used  by  the  military,  or  for  large  delay  multi-path  environ¬ 
ments,  true-time  delay  processing  is  needed.  The  length  of  the  delay  line  is  at  most 


reduced  below  the  noise  floor  of  the  system.  This  weight  distribution  in  the  complex 
plane  implies  that  each  of  the  received  signals  has  equal  importance.  For  adaptive 
algorithms,  the  adaptation  rate  is  inversely  proportional  to  the  number  of  weights 
and  the  computational  load  is  proportional  to  the  number  of  weights. 


This  thesis  presents  work  on  applying  lens  arrays  to  communication  and 
radar  systems.  It  covers  the  development  of  a  fast,  low-power  optically-controlled 
transmit/receive  lens  array  and  simulations  of  a  lens  array  adaptive  communication 
system.  Chapter  2  reviews  microwave  switch  design,  the  current  work  done  in  op¬ 
tical  switches  and  the  development  of  an  optically  controlled  resonant  microwave 


CHAPTER  2 


measure  of  how  well  the  switch  turns  “off”  the  load  and  is  expressed  as  l/|s2i|  for  impedance  in  the  “on”  path  and  a  very  large  impedance  in  its  “off”  path, 

the  switch  in  the  “off”  state.  Return  loss  appears  only  at  microwave  frequencies  and  As  mentioned  earlier,  in  microwave  circuits,  the  effective  impedance  of  a 

is  the  amount  of  power  reflected  at  the  switch  because  the  load  is  no  longer  ideally  device  depends  on  its  position  in  the  circuit.  The  shunt  configuration  for  SPDT 

matched  to  the  source.  It  is  expressed  as  l/|sn|  for  the  switch.  microwave  switch  exploits  this  property  to  improve  the  IL  or  ISO  of  the  switch 

in  some  cases.  For  the  “on”  path,  the  high  impedance  state  of  the  device  is  in 


an  electric  field  is  applied  between  the  metal  arm  and  the  pull  down  electrode  in  the  strength  and  the  stiffness  of  the  arm  (or  bridge),  the  arm  is  forced  to  the  contact 

gap.  The  electrostatic  forces  generated  pulls  the  arm  down  making  contact  with  quicker  and  relaxes  faster.  It  appears  that  the  MEMS  switches  are  fundamentally 

the  metal  contact  on  the  other  side  of  the  switch.  The  pull-down  electrode  usually  limited  to  a  “close”  time  of  about  1  ps  and  a  “open”  time  of  about  10  ps  due  to 

has  an  insulator  deposited  upon  it  to  prevent  the  shorting  of  the  control  voltage.  material  properties. 

By  reducing  the  electric  field  to  a  value  typically  less  than  the  threshold  voltage  Since  electrostatic  forces  are  used  to  control  the  MEMS  switch,  the  required 
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and  therefore  biases  a  microwave  PIN  diode.  The  third  technique  uses  a  optoelec- 
2.3  Developing  Photonic  Switch  tronic  device  to  control  the  bias  current  for  a  microwave  device,  Figure  2.3(c).  Each 

As  presented  in  the  previous  section  2.2,  most  of  the  work  done  in  photonic  technique  has  its  own  advantages  and  disadvantages,  and  has  been  investigate  to 

switches  has  been  done  in  order  to  develop  devices  that  could  sever  both  as  a  mi-  some  degree  in  this  project. 

crowave  and  optical  device.  In  this  work,  we  use  off-the-shelf  devices  for  both  the  The  first  technique  uses  a  optoelectronic  device  in  a  microwave  switch.  This 
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Figure  2.5:  Schematic  for  the  circuit  implementation  of  a  photo-MESFET 


Tkble  2.3.  The  performance  data  of  the  MESFET  and  two  photodiodes 
implement  the  photo-MESFET 


MicrostripTransmission  Line 


eral  switches  with  similar  enough  performance.  This  becomes  a  problem  when  the  P5tritro  0  rpi  ,  DT  ,  .  Tr  ,  x  ,  v 

Figure  2.14.  The  measured  RL  (a),  IL  (blue)  and  ISO  (red)  (b)  for  the  SPDT  i: 

resonant  frequency  of  the  switch  needs  to  be  the  same  as  the  resonant  frequency  of  Figure  2.13  before  (-)  and  after  (--)  the  new  diode  placement, 

the  patch  antenna  in  each  element  of  the  array. 


Microstrip  Line 


SPDT  Switch  Optical  Mount 


and  d=1.07  A/4  (yellow).  All  four  circuits  are  measured  with  the  same  calibration.  wittl  short  delivery  times,  and  require  only  a  single  supply  voltage.  Both  MMICs 

have  input  and  output  matching  filters  that  reject  the  lower  frequency  switching 
current,  Figure  2.23.  The  matching  filters  are  designed  assuming  gold  bond  wires 
that  are  0.254mm  (10 mil)  long  and  have  a  diameter  of  (0.7 mils).  The  MMICs  are 
mounted  on  a  brass  pedestal  to  minimize  total  bond  wire  lenght,  Figure  2.24. 
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Figure  2.31:  Illustration  of  microstrip  coupler  made  from  two  pieces  of  substrate. 


chip  PIN  diodes,  HPFs  and  LPFs. 


2.1  dB  ISO  for  the  series  configuration  and  0.3  dB  IL  and  2.3  dB  for  the  shunt  con-  1.5  nH.  The  IL,  ISO  and  RL  is  calculated  for  the  switches  paying  close  attention  to 

figuration,  Figure  3.3.  However,  the  T-network  configuration  performance  is  greatly  the  sensitivity  of  switch  performance  as  the  parameters  change.  The  lower  limit  for 

improved  through  reactance  cancellation  with  the  bond  wires  inductance.  the  inductance  at  1  nH  is  set  by  the  size  of  the  PIN  diode.  The  PIN  diode  sits  on 

There  are  three  design  parameters  that  affect  the  performance  of  a  PIN  a  0-5  mm  diameter  brass  rod  via,  limiting  the  gap  in  the  microstrip  line  to  1  mm. 


From  the  simulations,  the  IL  and  ISO  are  relatively  insensitive  to  changes  transmis¬ 
sion  line  impedance  and  length,  as  for  the  InH  case  in  Figures  3.4.  In  general, 
the  IL  and  ISO  improves  with  the  increase  of  transmission  line  impedance.  The  RL 
is  the  only  parameter  that  significantly  sensitive  to  design  parameters,  Figure  3.5. 
Overall,  1  nH  has  better  performance  than  the  1.5  nH,  however  caution  is  warranted 


Figure  3.12.  Measured  (blue  ~)  and  simulated  lumped  element  model  (red  --)  per¬ 
formance  of  HP F .  Measurements  above  5  GHz  are  calibrated  and  uncalibrated  below 
5  GHz. 


Bond  Wires 


Figure  3.17.  Three  different  test  configurations  used  to  measured  DC  voltages  and  Figure  3.19.  SPICE  simulation  of  Rise  time  and  fall  time  for  two  back-to-back  SPDT 

currents  for  developing  a  SPICE  model.  switches  used  in  the  unit  cell. 


Table  3.1.  Measured  input  and  output  power  for  back-to-back  SPDT  switch  showing 
maximum  input  power  for  the  SPDT  switch  is  at  least  17  dBm.  _ 
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the  setup  for  variations  in  cables  and  fiber  lengths. 

The  second  setup  is  for  low  optical  power  pulses,  which  requires  long  pulse 
widths,  Figure  3.27.  This  setup  uses  a  HP  8116X  function  generator  to  generate 
channel  1  a  square  wave  that  in  turn  drives  two  operational-amplifier  in  inverting  and  non- 


Optical  Power  (mW) 

Figure  3.34.  The  measured  rise  time  (red  •)  and  fall  time  (blue  •)  for  the  back- 
to-back  SPDT  switch.  The  measured  performance  is  consistent  with  the  simulated 
performance  (-)  and  within  the  simulated  error  bounds  (--). 


prototyping  milling  machine  mentioned  in  section  2.3.  The  milling  machine  can  mill  the  microwave  circuits. 

FR4  (PCB  material),  aluminum  and  brass,  as  well  as  the  softer  Teflon-based  sub-  4.1.1  Unit  Cell  Layout  The  layout  of  the  unit  cell  is  based  on  the 

strates.  Soft  metals,  such  as  copper  gum,  the  milling  tools,  while  low  temperature  layout  of  the  array,  such  that  the  performance  and  stability  of  the  unit  cell  is  as 

plastics  like  Nylon,  Delrin  and  polycarbonate  melt  from  the  high-speed  milling  tool  close  as  possible  to  that  of  the  array.  There  are  two  sides  to  the  unit  cell,  the  active 

friction.  component  side  and  passive  component  side. 


On  the  active  component  side,  effort  is  made  to  have  about  30  dB  of  isolation 
between  any  two  paths  that  can  carry  the  microwave  signal.  Since  the  maximum 
gain  of  the  MMICs  is  about  16  dB,  this  should  prevent  any  free  space  oscillations. 
The  electrical  and  mechanical  components  are  restricted  to  be  outside  an  area  of 


I 


The  lens  array,  Figure  4.8,  is  designed  as  a  cylindrical  lens  with  3  row  of  elliptical  in  shape,  it  is  impossible  to  get  a  perfectly  collimated  beam  using  a  single 

8  unit  cells.  The  slant  to  the  rows  allows  for  tighter  packing  of  the  unit  cells  with  lens.  It  is  important  to  collimate  the  laser  beam  since  the  optical  fibers  are  limited 

3/4  A  spacing  in  the  focusing  dimension  and  1 A  in  the  non-focusing  dimension.  The  to  an  acceptance  angle  of  25.4°  or  less,  determined  by  their  numerical  aperture  of 

focal  distance  to  diameter  ratio  is  1  and  the  points  of  perfect  focus,  section  5.1,  are  0.22.  The  optical  fiber  bundles  were  polished  by  hand  using  a  specially  made  jig 

at  ±45».  Figure  4.9  shows  the  simulated  radiation  patterns  for  the  lens  with  0»  and  to  hold  them  perpendicular  to  the  polishing  surface.  Unfortunately,  due  to  human 


Figure  4.10.  Simulated  normalize  radiation  patterns  for  the  0°  off  axis  feed  (red) 
and  the  30°  off  axis  feed  (blue)  in  the  focusing  dimension. 


l*4V  )  f(°-58J7  (  ,*3V)  (  l  4V)  (  MV  j  [  1.4V  ]  ( 1.4V 


4.4  Shaped  Pulse  Control  of  SPDT  Switch 

For  controlling  very  large  antenna  arrays,  the  optical  control  pulses  maybe 
further  optimized  in  shape  to  provide  more  efficient  use  of  the  optical  power.  As 
mentioned  in  section  3.1.4,  the  controlling  factor  on  how  fast  the  SPDT  switch 


unit  cells  to  be  controlled.  For  the  shape  optical  pulse  the  average  optical  power 
needed  per  unit  cell  is 
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where  the  /0  and  Iin  are  the  electric  current  density  for  the  antenna  and  input  current 
at  the  food  of  the  antenna  respectively. 

The  radiation  resistance  Is  essentially  a  normalization  of  the  total  power 


The  lens  object  contains  all  the  information  to  describe  the  lens:  size,  shape  (circular 
or  rectangular),  number  of  unit  cells  and  the  design  parameters  (F/D  and  Oq)  as  well 
as  a  complete  description  of  each  unit  cell  (gain,  delay,  feed  side  and  non-feed  side 
antenna  objects).  Since  the  lens  array  is  an  approximation  to  a  lens,  the  way  the 


Object  7:  Channel 

The  channel  describes  how  the  signals  for  a  user  is  received  by  the  array.  It  i 
collections  of  ncar-ficld  and  far-ficld  sources  that  represent  the  paths  that  the  us 
signal  took  to  the  array. 

Object  8:  QOLTF 


the  larger  the  number  of  unit  cells  in  the  array,  the  more  power  it  can  collect,  which 
argues  for  maximally  packed  unconstrained  lens  designs. 

5.3.2  Lens  Abberations  The  next  experiment  looks  at  aberration 
in  the  image  for  different  lens  designs  as  the  source  moves  off  axis.  For  an  ideal 
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Figure  5.17.  Image  patterns  for  sources  at  0°  (a),  20°  (b),  40°  (c)  and  60°  (d)  off 
axis  showing  the  aberrations  in  the  image. 
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Here  tho  £[•]  denotes  the  expectation.  Since  the  error  function,  <,  is  quadratic  in  Using  the  eigenvectors  of  R  to  rotate  the  system  to  its  principal-axis  coordinates 

W,  it  has  one  minimum  given  by  (where  the  signal  correlation  matrix,  R,  is  diagonal)  expressed  using  a  primed  vari- 

Wopt  —  R-1P  (A  7)  aklCj  ft  can  shown  that 

V'k  =  {I-2iiA)kV„'  (6.12) 

which  is  sometimes  known  as  the  Wiener  solution. 
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in  Equation  (6,12).  Then  the  noise  is  calculated  as 
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noise  to  the  output,  for  step  size,  p  =  . 

In  all  cases,  the  system  performance  improves  as  more  signals  are  added  to 
the  system,  which  is  expected.  The  added  signal  and  noise  powers  cause  the  weights1 
magnitude  to  adjust  inwards.  Figure  6.15  shows  the  change  in  the  optimal  weights 


inwards  to  minimize  the  amount  of  noise  at  the  output,  Equation  (6.47)  as  illustrated  <h“  incrftase<1  SNRx  means  a  decrease  in  total  noise  power  and  input  power  for  the 

by  Figure  6.14.  The  measured  performance  of  the  adaptive  system  is  slightly  less  8yStCm’  EqUtttion  (G'5)>  thc  gain  constant  is  lurgcr>  Equation  (6.14),  accounting  for 

that  that  of  thc  optimal  weight  system  showing  thc  effects  of  the  weight  noise  in  the  ***  incTease  ta  weight  n°ise-  The  varifttion  in  weight  noise  bfttween  the  ^fehts  is 

system.  Since  the  difference  between  the  optimal  weight  system  SNR  and  the  adap-  SmaU  consistcut  with  tllcory-  Thosc  simulations  still  raise  thc  question:  “Is  it 

tWo  weight  system  SNR  is  small,  the  weight  noise  does  net  contribute  considerable  P0SSible  *  ***  “  Signal  WhOSe  SNR  h  80  Sma11  *httt  the  addition  °f  this  signaJ  to 
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CHAPTER  7 


pattern  presented  in  section  6.2  appear  to  corresponds  with  the  beams  produce  by 
network  lenses.  Finally,  an  analytical  model  for  the  lens  array  transformation  should 
be  developed  that  describes  the  lens  array  loss,  aberrations  and  antenna  element  ef¬ 
fects  so  that  a  communication  model  can  be  developed. 
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The  Qi’s  shows  how  well  the  lens  coherently  focuses  the  re-radiated  signals  to  each  This  is  as  if  the  lens  array  does  not  focus  the  signals  but  illuminate  the  focal  surface 

detectors  and  are  constrained  such  that  equally.  In  this  case,  the  SNR  of  Equation  (7.28)  reduces  to 


Figure  7.9.  Image  pattern  for  “other”  user  at  -40°  (a),  0°  (b),  20°  (c)  and  25"  (d)  Figure  7.10.  The  movement  of  the  optimal  weight  for  lens  array  (a)  and  planar  array 

in  multi-user  simulation.  (b)  adaptive  systems  as  the  “other”  user  moves  closer  to  the  “desired”  user.  The 

final  realization  of  the  optimal  weights  arc  the  numbered  data  points. 


arrays  have  less  loss  than  lens  arrays,  These  results  suggest  that  lens  array  and  a 
planar  array  have  the  same  ability  to  null  out  “other”  users  and  that  there  is  not  a 
fundamental  limitation  in  achievable  radiation  patterns, 

7.3.4  Adaptive  Lens  Array  Systems  with  Multi-path  For  the 
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[20]  G.  Kriehn,  A.  Kirulata,  P.  E.  X.  Silveira,  S.  Weaver,  S.  Kraut,  K.  Wagner,  R.  T. 
Weverka,  and  L.  Griffiths,  “Optical  beamtap  beam-forming  and  jammer-nulling 
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function  plot. error (sub_non , window) 
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symbol_start=ones(lengtli(symbol)) ; 
symbol_start_opt=ones (length (symbol)) ; 
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function  [y_SNR] *y_SNR_ wiener (output , TF , desire.signal , N_s) 


